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Abstract 

We discuss some QCD constraints on light-cone tt meson wave 
function i/;(kj_,x) . The analysis is based on such general methods 
as dispersion relations, duality and PCAC. We calculate the asymp- 
totical behavior of the wave function (wf) at the end-point region 
(x — * 1 and — > oo) by analysing the corresponding large n— th 
moments in transverse (fcj_ n ) ~ n\ and longitudinal ((2a; — ~ 1/n 2 
directions. This information fixes the asymptotic behavior of wf at 
large k\ ( which is turned out to be Gaussian commonly used in the 
phenomenological analyses). 

We discuss one particular application of the obtained results. We 
calculate the nonleading "soft" contribution to the pion form factor 
at intermediate momentum transfer. We argue, that due to the spe- 
cific properties of ip(kj_,x), the corresponding contribution can tem- 
porarily simulate the leading twist behavior in the extent region of 
Q 2 : 3GeV 2 < Q 2 < 40GeV 2 , where Q 2 F(Q 2 ) ~ const. Such a mech- 
anism, if it is correct, would be an explanation of the phenomenologi- 
cal success of the dimensional counting rules at available, very modest 
energies for many different processes. We discuss some inclusive am- 
plitudes ( like Drell Yan and Deep Inelastic ) where intrinsic n meson 
structure might be essential. The relation to the valence quark model 
is also discussed. 
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1. Introduction 

The main goal of the present paper is the study of the hadronic wave 
functions with a minimal number of constituents within QCD. To be more 
specific, we are mainly interested in the k\ — behavior of the light cone wf 
in the transverse direction. 

The motivation for this interest is the following. As is known, at asymp- 
totically high energies the parametrically leading contributions to hard ex- 
clusive processes can be expressed in terms of the so-called distribution am- 
plitudes 4>(x), which itself can be expressed as an integral / ip(k^, x)d 2 k with 
nonperturbative wave function ip(kj_,x) , see review Jlp. Distribution ampli- 
tudes 4>{x) depend only on longitudinal variables X{ and not on transverse 
k\ ones. The same is true for inclusive reactions where structure functions 
depend on x, but not on k 2 ^. Thus, any dependence on k\ gives some power 
corrections to the leading terms. Naively one may expect that these cor- 
rections should be small enough already in the few GeV 2 region. However, 
as we argue later, this expectation does not seem work well in intermediate 
region. Thus, one can say that we study the pre -asymptotic behavior of the 
exclusive amplitudes. 

We shall find that ip(kj_, x) possesses the quite unusual properties, which 
lead to the broadening wf in the transverse directions. In terms of observable 
amplitudes it means that the characteristic scale is not ~ lGeV 2 (as naively 
one could expect), but lOGeV 2 . What is more important, this scale is not 
universal, but varies from process to process. 

The analysis of some inclusive reactions shows the same result- very of- 
ten the experimental data can not be explained within the standard scale- 
invariant description. An explicit introduction of some hadronic (phenomeno- 
logical) dimensional parameters into the structure functions is often required. 
Due to the fact that the structure functions can be expressed in terms of the 
same wave functions as ~ J \ip(k 2 L ,x)\ 2 d 2 k, we believe that the similar con- 
clusion (on importance of the pre-asymptotic behavior in the intermediate 
region) takes place for the inclusive amplitudes as well. 

As the simplest application of our wf we consider the ir meson form 
factor. Before to explain the qualitative results obtained in this paper, we 
would like to review a few important steps which have been taken in the 
investigation of the exclusive amplitudes. We hope that this short historical 
introduction will help to formulate the problem we want to address in the 
present analysis. 

3 The formal reason for that can be seen from the following arguments. At large 
energies the quark and antiquark are produced at small distances z ~ 1/Q — > 0, where Q 
is typical large momentum transfer. Thus one can neglect the z 2 dependence everywhere 
and one should concentrate on the one variable zQ ~ 1 which is order of one. One can 
convince oneself that the standard Bjorken variable x is nothing but Fourier conjugated 
to zQ. 
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In early seventies the famous dimensional counting rules were proposed 
0. The predictions of these rules agree well with the experimental data, 
such as the pion and nucleon form factors, large angle elastic scattering cross 
sections and so on. This agreement served as a stimulus for further theoretical 
investigations. The modern approach to exclusive processes was started in 
the late seventies and early eighties ||. We refer to the review papers 
for details . 

The main idea of the approach || is the separation of the large and small 
distance physics. At small distances we can use the standard perturbative 
expansion due to the asymptotic freedom and smallness of the coupling con- 
stant. All nontrivial, large distance physics is hidden into the nonperturba- 
tive wf in this approach. It can not be found by perturbative technique, but 
rather should be extracted from elsewhere. The most powerful analytical 
nonperturbative method for such problems is the QCD sum rules 0,0. 

The first application of QCD sum rules to the analysis of nonperturbative 
wf was considered more than decade ago M. Since then this subject is a 



very controversial issue [p|-||22| and we are not going to comment these quite 
opposite points in the present note. 

At the same time, the applicability of the approach at experimentally 



accessible momentum transfers was questioned []r0||,||13||. In these papers it 
was demonstrated, that the perturbative, asymptotically leading contribu- 
tion, is much smaller than the nonleading ("soft") contributions. Similar 
conclusion, supporting this result, came from the different side, from the 



QCD sum rules, pOflilPlHi where the direct calculation of the form factor has 
been presented at Q 2 < 3GeV 2 . This method, has been extended later for 
the larger Q 2 < lOGeV 2 [11|,[22]| with the same qualitative result: the soft 



contribution is more important in this intermediate region than the leading 
one. 

Now we are ready to formulate the question, which we want to address 
in the present paper. 

• If the asymptotically leading contribution can not provide the experi- 
mentally observable absolute values, than how can one explain the very good 
agreement between the experimental data and dimensional counting rules |, 
which are supposed to be valid only in the region where the leading terms 
dominate? 

It is clear, that the possible explanation can not be related to the specific 
amplitude, but instead, it should be connected, somehow, to the nonpertur- 
bative wave functions of the light hadrons (tt, p,p...) which enter the formulae 
for exclusive processes. The analysis of the tt meson form factor, presented 
below supports this idea. 



4 These rules unambiguously predict the dependence of amplitudes on dimensional 
parameters. In particular, Q 2 F 7T (Q 2 ) ~ constant., Q 4 F p (Q 2 ) ~ constant., s 7 ^-( , fp — > 
7r + n) = f(t/s) .. The experimental data are in a good agreement with these predictions 
in the large region of s, Q 2 at very modest energy and momentum transfer. 
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To anticipate the events we would like to formulate here the result of this 
analysis. The very unusual properties (which will be derived from QCD and 
not from quark model) of the transverse momentum distribution of the non- 
perturbative 7r meson wave function lead to the temporarily simulation 
of the dimensional counting rules by soft mechanism for the F^{Q 2 ) at the 
extent range of intermediate momentum transfer: 3GeV 2 < Q 2 < AOGeV 2 . 
In this region the soft contribution to F n (Q 2 )Q 2 does not fall-off, as naively 
one could expect, and we estimate it as F n (Q 2 )Q 2 ~ 0.3 -v- OAGeV 2 The 
leading twist contribution, after Sudakov suppression, gives, according to 
g,Pl,|2| a little bit less ( < 0.2GeV 2 ). 

• Therefore, our answer on the formulated question is the following. The 
nonperturbative wave functions possess (along with the standard small scale 
k\ ~ (330MeV) 2 ) the new, larger scale (~ lGeV 2 ). Precisely this new scale 
defines the regime where the asymptotical formulae start to work. 

We believe that the same features of the wf may affect the analysis 
of inclusive amplitudes as well where some pre-asymptotic effects might be 
essential. 

The paper is organized as follows. In the next section we define the 
nonperturbative wf through its moments. We focus on the properties of 
the two particle leading twist wf and its quark longitudinal and transverse 
distributions. We recall our previous analysis regarding the nonperturbative 
wf and formulate the main constraints which have been obtained from QCD 
analysis. In section 3 we model the wf which satisfies these constraints. 

Section 4 is devoted to the calculation of the soft contribution to pion 
form factor based on the model wfs obtained in the previous section. Let 
us stress from the very beginning: we are not pretending to have made a 
reliable calculation of the form factor here. We discuss some very general 
properties of the amplitudes which are related to the specific features of the 
ip{k 2 L , x). We illustrate how these features change the behavior of the form 
factor in intermediate region of Q 2 . Section 5 is our conclusion and outlook. 

2. Constraints on the nonperturbative wave function %p(k 2 L) x). 

First of all let us review some essential definitions and results about non- 
perturbative wf. We define the pion axial wave function in the following 
gauge-invariant way: 

iM A {zq,z 2 ) = (0\d(z) W e l9 ^ A ^u(-z)\n(q)) (1) 
= £ -<0|d(0) 7M 75(«* A,)M0)|7r(g)>, 

n n - 

where D V =D V — D v and iD^ = id^ + gA a ^ is the covariant derivative. 
From its definition is clear that the set of different tt meson matrix elements 
defines the nonperturbative wave function. 



3 



The most important part (at asymptotically high q 2 ) is the one related 
to the longitudinal distribution. In this case z 2 ~ the wf depends only 
on one zq- variable. The corresponding Fourier transformed wave function 
will be denoted as 0(£) and its n— th moment is given by the following local 
matrix element: 

(0\d lvl5 (i D, zXu\v(q)) = iUq v {zq) n (C) = iUq v {zq) n (2) 



— q 2 — > oo, zq ~ 1 £ = x\ — X2, x\ + £2 = 1, z 2 = 0. 

Therefore, if we knew all matrix elements (@) ( which are well-defined ) we 
could restore the whole distribution amplitude </>(£)• The QCD sum rules ap- 
proach allows one to find the magnitudes only the few first moments ||. As is 
known, this information is not enough to reconstruct the wf; the parametric 
behavior at £ — > ±1 is the crucial issue in this reconstruction. 

To extract the corresponding information, we use the following duality ar- 
gument. Instead of consideration of the pion wf itself, we study the following 
correlation function with pion quantum numbers: 

i ( dxe^ x (0\TJl(x), Jo(0)|0) = (zq) n+2 I n (q 2 ), jj| = d lv z vl5 {i D, ^) n «(3) 



and calculate its asymptotic behavior at large q 2 . The result can be presented 
in the form of the dispersion integral, whose spectral density is determined 
by the pure perturbative one-loop diagram: 



ImIn{sr rt = 3 

tt Jo s — q 2 4n(n + l)(n + 3) 

We assume that the tt meson gives a nonzero contribution to the dispersion 
integral for arbitrary n and, in particular, for n — > oo. Formally, it can be 
written in the following way 



1 r s " 1 r°° 

- / dslml(sy n ert = - / dslml(sr n , (5) 

TT JO TT JO 



Our assumption means that there are no special cancelations and tt meson 
contribution to the dispersion integral is not zero, i.e. £"(11) ^ 0, where 
we specified the notation for the longitudinal distribution. In this case at 
q 2 — ► oo our assumption (H) leads to the following relation: 



A 2 e>(™-°o)-^r ( 6 ) 



3S. 
47T 2 n 

It unambiguously implies the following behavior at the end-point region 



(o = £ dtem - i/n 2 , - ±i) - (i - a- (7) 
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Few comments are in order. We consider the nonperturbative correlation 
function (^). Thus the behavior (g) should be fulfilled for any nonperturba- 
tive wf. The perturbative as well as nonperturbative corrections will change 
the duality interval S%(\\) in the formula M) in comparision with perturbative 
one-loop calculation. However, l/n 2 - behavior remains unaffected. 
Thus, our first constraint looks as follows: 

.1 0(£_>±1)_,(1_£2)_ 

We want to emphasize that we did not use any numerical approximation 
in this derivation. Therefore, the constraint has very general origin and 
it should be considered as a direct consequence of QCD. Only dispersion 
relations, duality and very plausible assumption formulated above have been 
used in the derivation (»1). 

Now we want to repeat these arguments for the analysis of the transverse 
distribution. To do so, let us define the mean values of the transverse quark 
distribution by the following matrix elements: 

-> (2n - IV - 

(0|d 7 „7fe(i^ g to t.K?)> = iUqu(-t 2 ) n { { 2n)\\ {kr) - (8) 

where D v is the covariant derivative, acting on the one quark and trans- 
verse vector = (0, t, 0) is perpendicular to the hadron momentum = 
(<?Cb 0j_, q z ). The factor ^T^jj; is introduced to (H) to take into account 
the integration over angle in the transverse plane: / d<p(cos <fi) 2n / f d<p = 
(2n-l)!!/(2n)H. 

We interpret the (k\) in this equation as a mean value of the quark 
perpendicular momentum. Of course it is different from the naive, gauge 
dependent definition like (0| d r y u 'y 5 d'l_u\ii(q)) , because the physical transverse 
gluon is participant of this definition. However, the expression (§) is the 
only possible way to define the k\ in the gauge theory like QCD. We believe 
that such definition is the useful generalization of the transverse momentum 
conception for the interactive quark system. 

To find the behavior (fc^ 1 ) at large n we can repeat our previous duality 
arguments with the following result fj 

Z, 2 ^)^ ~ n\ => ~ nl (9) 

This behavior has been obtained in ref. ||24|| by analysing the perturbative 
series of the specific correlation function at large order. The dispersion rela- 
tions and duality arguments translate this information into the formula (H). 
Let us repeat again: any nonperturbative wave function should respect eq.(|9|) 
in spite of the fact that apparently we calculate only the perturbative part ( 
see comment after formula (0)). 

5 Here and in what follows we ignore any mild (nonfactorial) n-dependence. 
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The nice feature of is its finiteness for arbitrary n. It means that the 
higher moments 



do exist. In this formula we introduced the nonperturbative wf ip(kj_,£), 
normalized to one. Its moments are determined by the local matrix elements 
(H). The relations to Brodsky and Lepage notations ^ bl(xi, k±) QXJ and to 
longitudinal distribution amplitude 0(£) introduced earlier, look as follow: 

f 167r 2 r -> r 1 

^BL(x u k ± ) = ^-^^,k ± ), j dkiij(ki,o = m, J_ x dm) = xio) 

where f„ = 133MeV. The existence of the arbitrary high moments (fcj_ n ) 
means that the nonperturbative wf, defined above, falls off at large trans- 
verse momentum k\ faster than any power function. The relation (|||) fixes 
the asymptotic behavior of wf at large k\. Thus, we arrive to the following 
constraint: 

•2 (kf) = J dk 2 ± d^k 2 ^(k 2 ± , ~ n\ n -> oo. 

We can repeat our duality arguments again for an arbitrary number of 
transverse derivatives and large (n — > oo) number of longitudinal derivatives 
with the following result f23|j : 

.3 Jdklk 2 Mkl^-±l)~(l~e) k+1 - 
For the k = we reproduce our previous formula for the function: 
(j>(£ -> ±1) = fdk 2 ± ip(kl^ ±1) ~ (1 - £ 2 ). The constraint («3) is 
extremely important and implies that the k\ dependence of the fp{k\,^) 
comes exclusively in the combination k\j{\ — £ 2 ) at £ — > ±1. The 
byproduct of this constraint can be formulated as follows. The standard 
assumption on factorizability of the ip(kj_,^) = ^(^1)0(0 does contradict 
to the very general properties of the theory. Thus, the asymptotic behavior 
of the wf turns out to be Gaussian one with the very specific argument: 

^(kl - oo,0 ~ exp(--M_.), (ii) 

We would like to pause here in order to make the following conjecture. The 
Gaussian wf (reconstructed above from the QCD analysis) not accidentally 
coincides with the harmonic oscillator wf from constituent quark model. To 
make this conjecture more clear, let us recall few results from the constituent 
quark model. 



It is well known [55] that the equal- time wave functions 

?PcM(q 2 ) ~exp(-g 2 ) (12) 

of the harmonic oscillator in the rest frame give a very reasonable description 
of static meson properties. Together with Brodsky-Huang-Lepage prescrip- 



tion [26 1 , [27 connecting the equal -time and the light-cone wave functions 
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of two constituents (with mass m ~ 300Me\^) by identification 



-2 k± + 777,2 2 ; / -2\ / / ^1 + 777,2 2\ 

one can reproduce the Gaussian behavior ( |TH ) found from QCD. It means, 
first of all, that our identification of the moments (§) defined in QCD with 
the ones defined in quark model, is the reasonable conjecture Q. 

However, there is a difference. In quark model we do have a parameter 
which describes the mass of constituent m ~ 300MeV. We have nothing like 
that in QCD. This difference has very important consequences which will be 
discussed later. 

We would like to put a few more constraints on the list. But before to 
do so, we have to emphasize the difference between constraints («1 — «3) 
discussed above and the ones which follow. The first three constraints have 
very general origin. No numerical approximations have been made in the 
derivation of the corresponding formulae. The only what have been used 
are dispersion relations, duality and very plausible assumption formulated 
above. The main idea of the derivation of all these constraints is one and the 
same: we calculate some correlation function in QCD using the asymptotic 
freedom. The dispersion relations and duality arguments transform these 
properties into the constraints on hadronic matrix elements. 

The constraints we are going to discuss now have absolutely different 
status. They are based on the QCD sum rules with their inevitable numerical 
assumptions about higher excited states in QCD. Thus, they must be treated 
as an approximate ones. The well-known constraint of such a kind is the 
second moment of the distribution amplitude in the longitudinal direction 

& 

.4 (a = / dtm? - OA 

(The asymptotic wf corresponds to (£ 2 ) = 0.2 ). Such a result was 
the reason to suggest the "two-hump" shape wf || which meets the above 
requirement. The number cited as the constraint («4) has been seriously 



criticized in refs. P|-|f[2]. The point for criticism was exactly the assumption 
about the role of the excited states in the sum rules. We can not answer 
on this criticism within standard QCD sum rules approach. Thus, in what 
follows we shall discuss both possibilities: the narrow (asymptotic) wf and 
the wider one (with larger (£ 2 ) > 0.2). 

The next "numerical" constraint is the second moment of the wf in the 
transverse direction defined by equation @ and calculated for the first time 

6 The same method can be applied for the analysis of the asymptotical behavior of the 
nucleon wf which in obvious notations takes the form: 

^nucleonik J_i * °°7 x i) ~ ex P( — / y ") 

3jA 
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in || and independently (with quite different technique) in [p§| . Both results 
are in a full agreement to each other: 

•5 (fcl) = i {J ^n§^ - ^ - O.lGeV 2 , m 2 ~0.8GeV 2 . 

Essentially, the constraint («5) defines the general scale of all nonpertur- 
bative phenomena for the pion. It is not accidentally coincides with 300MeV 
which is the typical magnitude in the hadronic physics. 

All numerical values obtained within QCD sum rules approach correspond 
to the normalization point ~ lGeV 2 . To model the wf we need to know their 
renormalization properties. The anomalous dimensions of the longitudinal 
operators are well known, see e.g. [J], ||. For our particular case we can write 
it in the following way: 

The anomalous dimensions of the higher twist operators (which are re- 
lated to transverse moments) are less familiar. It has been calculated for the 
operator describing the mean value of the transverse momentum k\ \ll 



(kiU = • (kiU. (i4) 

To study the fine properties of the transverse distribution it is desired to 
know the next moment. The problem can be reduced to the analysis of the 
mixed vacuum condensates of dimension seven |[23|| : 



/giv _ 1 [ -3(qg a^G^ax^G^q) 13(gg G^G^q) 

[ ±} 8 1 4{qq) 9(qq) 1 1 

We analyzed the magnitudes for these vacuum condensates with the following 
result: the standard factorization hypothesis does not work in this case. The 
factor of nonfactorizability K ~ 3.0 4- 3.5 [25], [29|. The eq.([r5|) defines the 



new numerical constraint on the transverse distribution. We prefer to express 
this constraint not in terms of the absolute values, but rather, in terms of 
the dimensionless parameter R which is defined in the following way: 

•6 R= ~ 3K ■ <g2G ^ 4 G ^ } -5 4-7, ml ~ 0.8GeV 2 , 

<fei> 2 m 

where we use the standard values for parameter wig and gluon conden- 
sate 0. We would like to emphasize that the fluctuations of the transverse 
momentum are large enough. The quantitative characteristic of these fluc- 
tuations is parameter R^> 1. In terms of wave function this property means 
a very unhomogeneous distribution in transverse direction. 

As we already mentioned, to model wf we need to know the renormaliza- 
tion properties of the higher twist operators. Unfortunately, we do not know 
them (except for the lowest one, (0)). However, one can argue, that in the 
large N c limit, the main contribution to anomalous dimensions can be found 
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from the formulae like (0), where matrix elements are expressed in terms of 
vacuum condensates. In the same iV c — > oo limit the renormalization prop- 
erties of these condensates (but not their absolute values) can be estimated 
(with the same accuracy) by applying the factorization procedure. Thus, 
one could expect that this prescription gives a reasonable numerical accu- 
racy (~ 1/N C ) for the renormalization properties of higher twist operators. 
However, the exact calculations are highly desired and welcome ||30|| . 

We would like to check this prescription for the operator with known 
dimension (|14j). The anomalous dimensions for the chiral condensate and 



mixed condensate are known: 



'oc s (li 2 ) 



{qqU = ■ m„ (16) 



Thus, our relation («5) gives the following prescription for the evolution 
formula under the renormalization group transformation: 



(qiga^G 1 \ u ^q) a s {fn) 30 
(qq) a 8 (// 2 )' 



(kl),, ~ ^J^AIL „ (^)ff ■ (kl),,, (17) 



instead of exact formula (|1J). We consider this good numerical agreement 
as a justification for the analogous estimation for different operators with 
unknown anomalous dimensions. In particular, we expect that the matrix 
element (0) is not changed strongly under the renormalization group 

transformations because after the factorization this operator reduces to the 
gluon condensate which is renormalization invariant. 

At the same time our dimensionless parameter R is changed strongly and 
we estimate it as follows: 

Let us summarize the results of the section. The constraints («1 — «3) 
have very general origin and should be fulfilled in any phenomenological 
description based on QCD. The numerical constraints («4 — «6) have much 
less generality because they have been obtained from QCD sum rules with 
inevitable for this method approximations. All numerical results obtained 
from QCD sum rules are normalized at fi 2 ~ lGeV 2 . At the same time, 
the model wf we are going to construct should be normalized at the lowest 
possible point which is about fi 2 ~ 0.25GeV 2 . We shall use the anomalous 
dimensions shown above in order to evaluate all numerical constraints at the 
lowest normalization point. 
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3. The model wave function. 



Let us start our discussion from the analysis of the wf motivated by 
constituent quark model ]TI|, |2"7[ | (CQM) []. Such a function is known 
to give a reasonable description of static hadron properties. The Brodsky- 
Huang-Lepage prescription |27| leads to the following form for the pion wf: 

^,. W = ^exp(-J|±^_), (19) 

We call this function as the constituent quark model wf. It satisfies two 
constraints («2, «3), but not (»1) because of the nonzero magnitude for the 
constituent mass m. We take the standard set for QCD parameters: 

a s (lGeV) = 0.34, A QCD = 200MeV. (20) 

The lowest possible normalization point /io is defined as the place where 
a s (fj, ) ~ 0.7. This corresponds to the QCD sum rules analysis 0, where 
"almost" renormalization invariant combination a s (qq) 2 comes into the game 
and it is numerically well known. At the same time the chiral condensate 
(qq) at the lowest possible normalization point is known from PCAC. 

We made the standard choice for the parameter m ~ 330MeV in accor- 
dance with its physical meaning. The parameter (3 is determined from the 
numerical constraint («4) for the mean value (k\) f\. Parameter A is deter- 
mined by the normalization eq.flTUp. As we mentioned earlier, we have to 
renormalize all moments to the lowest possible point to model wf. In our 
case the evaluation of (k\) is defined by eq. (|I~4"D . With our set of parameters 
( p0|) and («5) we have the following mean value for h\ at lowest normalization 
point: 

(k 2 ± )^ = O.UGeV 2 , (21) 

which will be used through this paper. This number corresponds to f3 ~ 
0.3GeV which is within reasonable parametric region. To make function 
wider in the longitudinal direction (constraint «4) one can insert to formula 

7 Here we neglect all terms in QCM related to spin part of constituents. In particular, 
we do not consider Melosh transformation and other ingredients of the light cone<=>equal 
time connection. It does not effect qualitative results presented in the next section. 

8 Do not confuse our parameter (kj_) which is well defined in terms of QCD matrix 
element (||) with the one which is defined in terms of CQM as follows: 

,p > _ = / d 2 k 2 L dx\ip(k 2 ± , x)cQNi\ 2 k]_ 
J d 2 k\dx\'ip{k\ 1 x)cQM\ 2 

and should be extracted from somewhere else. Numerically they are not very different, but 
we prefer to use one and the same procedure to specify parameters for all wave functions. 
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([Op the additional factor 



(1 + gM[{2x - l) 2 - 1]), ffOii) = • (22) 

with additional parameter g(fi). With this new parameter <?(//) one can 
adjust (£ 2 ) as appropriate. For the asymptotic distribution amplitude pa- 
rameter g — 0. 

We are ready now to discuss the model wf in QCD. Before to design 
"0qcd(^1,^), let us explain what do we mean by that. We define the non- 
perturbative wave function ip(k±, x, h)qcd through its moments which can be 
expressed in terms of the nonperturbative matrix elements ([l|). As is known, 
all nonperturbative matrix elements are defined in such a way that all gluon's 
and quark's virtualities smaller than some parameter \i (point of normaliza- 
tion) are hidden in the definition of the " nonperturbative matrix elements" . 
All virtualities larger than that should be take into account perturbatively. 
In particular, all perturbative tails like l/kj_ should be subtracted from the 
definition of the nonperturbative wf. The same procedure should be applied 
for the calculation of nonperturbative vacuum condensate (G 2 ^), where the 
perturbative part related to free gluon propagator l/k 2 should be subtracted. 

With these general remarks in mind we propose the following form for the 
nonperturbative wave function if)(lc±, x, Ho)qcd at the lowest normalization 
point: 

^(kl,x,fi ) Q cD = Aex P(~ ^7fZ^y ) ' i 1 + 9(M(2x - l) 2 - |]}.(23) 

In comparison with the constituent quark model the "only" difference is the 
absence on the mass term ~ m in the exponent. As we shall see in the next 
section it does make a big difference. This function satisfies all fundamental 
constraints («1 — «3). The dimensional parameters can be determined from 
the numerical relations («4 — «5) in the same way as before. Parameter (3 for 
this parametrization is found to be (3 ~ 0.3 GeV (it corresponds to R = 2.2 
and (k 2 ± ) = O.UGeV 2 ) 

We would like to emphasize that the nonzero mass in ip(kj_, x)cqm was 
unavoidable part of the wave function. We do not see any room for such 
term in QCD, because its presence would mean the following behavior of the 
large moments in longitudinal direction: 

(O = J\CH0 ~ /_ 1 i ^Cexp(- T 4^) ~ exp(-v^), n - oo. (24) 

It is in contradiction to 1/n 2 behavior (|7p, («1). Let us stress that such a 
behavior is the result of the calculation of the correlation function (|3|). We 
do not see any possibilities to change this behavior from 1/n 2 to exp(— \/n). 
If such a behavior were occurred, it would mean that the strong cancellation 
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between one-loop diagram (f|) and higher loop corrections ~ a k s takes place. 
Such a cancellation looks even less probable, if one takes into account that 
the aforementioned cancellation must take place for each given number n at 
large n. We do not believe that it might happen in QCD. 

Up to now we did not discuss the influence of our last "numerical" con- 
straint («6) denoted by R. Large number for this parameter means a notice- 
able fluctuations of the momentum in transverse direction. To satisfy this 
constraint we need to spread out the wf to make it wider. It can be done in 
arbitrary way. In particular, one may try to put one more hump apart from 
the main Gaussian term described by eq . (|23|) . The only requirement is: it 
has to fall off fast enough at large k\. 

With these remarks in mind we suggest the following QCD motivated 
wf (we call it as ip(kj_, x, Ho)qcd+) which can be adjusted to satisfy all six 
requirements mentioned in the previous section: 

i[)(k 2 x , x, ^)qcd+ = A{e~*P*ti-*) + c • e^s^u-*) -02 }. (25) 

{l + g^Wx-iy- 1 -}}. 

The physical meaning of the parameters c, / is clear. Parameter c determines 
the magnitude of the second hump and parameter / describes its distance 
from the main term. As we shall see in order to match parameters I, c with 
the calculated ratio R we need to have the magnitude of the second hump 
about 1/10. The maximum of the second hump is located in k\ ~ 0.8Gey 2 
region. To be more specific, we renormalized the parameter R found from 
QCD to the lowest normalization point according to formula flT8|). It is found 
to be 

i?(/V)-3-4 (26) 

We display this function with parameters (3 = 0.15 GeV , c = 0.15, 1 = 30 
(which correspond to R ~ 4, (k\) = O.lAGeV 2 ) on Fig.l at the central point 
x = 1/2. 

Let us summarize. We constructed three wave functions. The first one, 
ipCQM is motivated by quark model with its specific mass parameters. Two 
other models are motivated by QCD consideration. Despite of this difference, 
all these models have Gaussian behavior at large k\. However, in the case of 
ipCQM this behavior is related to the nonrelativistic ocsilator model, while for 
QCD motivated models this behavior is provided by constraints discussed in 
the previous section. 

Contrary to the CQM, the QCD motivated wave functions do not contain 
the mass parameter m ~ 300MeV which is an essential ingredient of any 
quark model. Such a term is absolutely forbidden from the QCD point of 
view. 
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The difference between i^qcd and 4>qcd+ is not fundamental, but rather 
quantitative. Nevertheless, we believe that it is worth to mention some new 
effects (noticeable fluctuations of transverse momentum) which the function 
i/jqcd+ brings. The broadening in the transverse direction is the main differ- 
ence between functions ^qcd and iPqcd+- In the next section we discuss the 
contribution to the pion form factor caused by these three wave functions. 
We shall see the qualitative difference in behavior on Q 2 , which is our main 
point. 

4. Pion Form Factor. 

The starting point is the famous Drell-Yan formula [[31] (for modern, 
QCD- motivated employing of this formula, see |[J), where the F^Q 2 ) is 
expressed in terms of full wave functions: 

Ftt{Q 2 ) = I ^^HT BL {x, k ± + (1 - x)q ± )* BL (x, k ± ), (27) 

where q 2 = —q±2 = —Q 2 is the momentum transfer. In this formula, the 
^ bl(x, k±) is the full wave function; the perturbative tail of ^ bl(x, k±) 
behaves as a s jk\ for large k\ and should be taken into account explicitly in 
the calculations. This gives the one-gluon-exchange (asymptotically leading) 
formula for the form factor in terms of distribution amplitude <p(x) ||. 

Below is the QCD motivated interpretation of this formula. Let us re- 
mind, that the formula ( p7|) takes into account only the valence Fock states. 
The formula would be exact if all Fock states were taken into account. Be- 
sides that, k\ in this formula was originally thought to be the usual (not 
covariant) perpendicular momentum of the constituents, and not the mean 
value (k 2 ^) defined in QCD, as a gauge invariant object. However we make 
the assumption that it is one and the same variable. The physics behind 
of it can be explained in the following way. 

In the formula (|27|) we effectively take into account some gluons (not all 
of them), which inevitably are participants of our definition of wf. These 
gluons mainly carry the transverse momentum (which anyhow, does not ex- 
ceed QCD scale of order ~ lGeV) and/or small amount of the longitudinal 
momentum. The contributions of the gluons carrying the finite longitudinal 
momentum fraction are neglected in (|27|) . This is the main assumption. It 
can be justified by the direct calculation Q of quark-antiquark-gluon (with 
finite momentum fraction) contribution to tt meson form factor at large Q 2 
within the standard technique of the operator product expansion. By techni- 
cal reasons the corresponding calculation has not been completed, however it 
was found that the characteristic scale which enters into the game is of order 
lGeV 2 . Thus, it is very unlikely that these contributions can be important at 
Q 2 3> lGeV 2 . The second calculation, which confirms this point, comes from 



the light cone QCD sum rules p2|. This is almost model independent calcu- 
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lation demonstrates that the quark-antiquark-gluon (with finite momentum 
fraction) contribution does not exceed 20% at available Q 2 . 

Thus, we expect, that by taking into account the only "soft" gluon contri- 
bution (hidden in the definition of k\ (Q)), we catch the main effect. Again, 
there is no proof for that within QCD, and the only argumentation which 
can be delivered now in favor of it, is based on the intuitive picture of quark 
model, where current quark and soft gluons form a constituent quark with 
original quantum numbers. No evidence of the gluon playing the role of a 
valent participant with a finite amount of momentum, is found. 

From the viewpoint of the operator product expansion, the assumption 
formulated above, corresponds to the summing up a subset of higher- 
dimension power corrections. This subset actually is formed from the in- 
finite number of soft gluons and unambiguously singled out by the definition 
of nonperturbative wf (0). 

In the following, we preserve the notation ^bl{x, k±) for the nonpertur- 
bative, soft part only. It should not confuse the reader. 

The formula Q2"7|) is written in terms of Brodsky and Lepage notations 
Jl[]; the relation to our wave function ip(k±, x) is given by formula ([10|). 

With these general remarks in mind, we would like to present the results of 
calculation, based on three wave functions discussed in the previous section. 

The first calculation, based on ipcQM(k±,x) is the standard one. The 
analogous calculations with ocsilator-like wf have been done many times 
with many additional improvements, see i.e. ]nj. One can fit the dimensional 
parameters in such a way that the description at low Q 2 will be perfect. 
However, our goal is different and we are interested in the behavior at large 
enough Q 2 ^> lGeV 2 . We expect that described here approach makes sense 
only at high enough Q 2 . We display the corresponding behavior as a curve 
1 on Fig.2. 

The main feature of this behavior - it gives very reasonable magnitude 
for the intermediate region about few GeV 2 and it starts to fall off right after 
that. We expect that any reasonable, well localized, based on quark model 
wave function with the scale ~ (k±) ~ m 2 leads to the similar behavior. Let 
me stress: we are not pretending to have made a reliable calculation of the 
form factor here; we displayed this contribution only for the illustrative 
purposes. 

Currently, much more interesting for us is the calculation, based on 
QCD motivated models. We display the corresponding contribution to the 
Q 2 F n (Q 2 ) for the i{)qcd{x, k±) on Fig. 2 as curve 2. The qualitative dif- 
ference between this curve and the previous one (curve 1), is much slower 
fall off at large Q 2 for the model wf i^qcd- The qualitative reason for that 
is the absence of the mass term in i^qcd, see discussion after the formula 
(p4|). Precisely this term was responsible for the very steep behavior in the 
previous calculation with quark model wave function. 

The declining of the form factor getting even slower if one takes into 
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account the property of the broadening of wf in transverse direction. The 
corresponding contribution based on iPqcd+ is displayed on Fig. 2 as curve 
3. We notice an additional slowing down of the declining of the magnitude 
Q 2 F(Q 2 ) in the intermediate region of Q 2 . The explanation of this effect 
is the following. When Q 2 is getting bigger and bigger, the contribution 
coming from the overlap of two humps (i/jqcd+ at x = 1/2 is displayed on 
Fig.l) starts to grow. These humps are well separated (in order to satisfy 
constraint (»6)) from each other by the value of order lGeV 2 . Therefore, we 
expect that this contribution starts to grow at high enough Q 2 3> iGeV 2 . 
As the result, the curve 3 looks more horizontal than the previous one. It 
is getting lower because the general scale has been changed when we passed 
from iPqcd to ipQCD+- Let us remind that we keep (k±) fixed (»5) in all 
cases. It leads to some changes of dimensional parameters because the small 
hump in the ipQCD+ gives a noticeable contribution to (k±) in spite of the 
fact that it comes with very small relative weight (coefficient c ~ 1/10). 

Our last qualitative remark is some note that the results strongly depend 
on parameter (£ 2 ). We displayed on Fig.3 the same three curves as on Fig.2 



with the only difference in coefficient g(fj,) (|22|, [23|, |25|) We set g([J>o) = 2 which 
corresponds to (£ 2 ) = 0.3. It makes a wf wider in longitudinal direction. The 
soft contribution getting bigger when a wider (in longitudinal direction) wave 
function is used. The same effect was observed in the recent calculation |22|] . 
where a quite different method has been used. 

The contribution under consideration is subject to Sudakov corrections. 
An estimate of these corrections reveals that they are small enough in this 
intermediate Q 2 region. Besides that, we will be on the safe side if we say 
that the hard (leading twist ) contribution to Q 2 F(Q 2 ) ~ 0.2GeV 2 pfl, PH, 



22]] . It should be added to the soft terms displayed on Fig.2, Fig.3. 

The precise fitting of the pion form factor was not among the goals of this 
paper. Rather, we wanted to demonstrate how the qualitative properties of 
a nonperturbative wf, derived from the QCD analysis might significantly 
change its behavior. 

5. Summary and Outlook. 

The main goal of the present paper was the analysis of the nonperturba- 
tive wf from QCD point of view. We found qualitatively different results in 
comparison with the wave functions motivated by quark model. We believe 
that this difference is responsible for the qualitative explanation of dimen- 
sional counting rules which work well even at very modest energies. 

The standard point of view for the phenomenological success of the di- 
mensional counting rules is the predigest that the leading twist contribution 
plays the main role in most cases. We suggest here some different explana- 
tion for this phenomenological success. Our explanation of the slow falling 
off of the soft contribution with energy is due to the specific properties of 
nonperturbative wf. In particular, we found a new scale (~ lGeV 2 ) in the 
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problem, in addition to the standard low energy parameter (k±) ~ O.lGeV 2 . 

Our next remark can be formulated as follows. Exclusive, as well as in- 
clusive amplitudes can be expressed in terms of the one and the same 
particular hadron wf. Therefore, if our explanation (related to specific form 
of i])(k\,x)) of a temporary simulation of the leading twist behavior is con- 
sidered reasonable one , then: 

1. in the analysis of inclusive amplitudes one may expect the same effect 
( it is our conjecture ); 

2. one may try to implement the intrinsic transverse momentum depen- 
dence into the inclusive calculations. 

In particular, one may try to use the following prescription for the 7r 
meson distribution function (and anologously for nucleon) at x =>■ 1: 

G q/ «(x, Q 2 )^{J ex P(-^(TT^)} G ^' ^ ( 28 ) 

The analogous formula (without x dependence in the exponent) has been 
suggested many years ago [p2]| , see also [ 33f| . Our remark is that this x 



dependence is essential point and should be introduced to the formula to 
satisfy our constraints. In terms of H32] , [ 33|| it corresponds to the non- 



universality of their "constant" (k±) which now will depend on x. 

To support this conjecture, we would like to mention few inclusive pro- 
cesses where the intrinsic transverse distribution might be essential. First of 
all, it is Drell-Yan amplitude tc + iV — > fi~[i + + X which is parametrized as 



follows (for references and recent development see |3~| ) : 

~ 1 + ^ cos 2 9 + fj, sin 29 cos <fr + - sin 2 9 cos 20. (29) 
a ail 2 

Here 9, <p are angles defined in the muon pair rest frame and A, /i, v are 
coefficients. In the naive parton model the coefficient are A = = v = 
0. Experimental results do not support this naive prediction. Recently, 
some improvements have been made |53], but some problems are still remain. 



In particular, the Lam- Tung sum rule p5| , 1 — A — 2v = is violated by 
experimental data and the improved model |33| still can not explain the 
behavior 1 — A — 2u as a function of Q\ ( Q 2 ± is the transverse momentum 
of the lepton pair). 

Due to the fact that Q 2 is not large enough in this experiment one may 
expect that the intrinsic distribution (|2"%| ) might be essential. 

Analogously, we would expect that the n meson structure (p8|) may effect 
the analysis of azimuthal asymmetries in semi-exclusive amplitudes like I + 
p — > I'+h+X. For references and recent development see f56fl , where intrinsic 
k\ has been introduced in the standard way without x dependence in the 
exponent. 

One may find many examples like that where the standard parton picture 
does not work well. We would like to mention here the recent analysis [^7 
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of the direct photon production (n + p —>■ 7 + X), with the result that per- 
turbative QCD can not explain the data. Some nonperturbative broadening 
factor in transverse direction should be implemented. One may hope that 
formula may improve the agreement with experiment. 
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FIGURE CAPTIONS 



Fig.l QCD Wavefunction ( unnormalized ) ip(kj_, x, hq)qcd+ versus trans- 
verse momentum k±. 

Fie. 2 Pion Form Factor Q 2 ■ F W (Q 2 ) versus Q 2 . Line 1 corresponds to the 
il)(kl,x) C QM, with parametres, R ~ 2.0 ©, (k 2 ± ) = 0.14 GeV 2 (|) , line 
2 corresponds to the i/j(kj_, x, Hq)qcd with parametres, R ~ 2.2, (kj_) =0.14 
GeV 2 , (C, 2 ) =0.2, and line 3 shows the result for two humped ( in transverse 
direction ) wavefunction i/j(kj_, x, /io)qod+, with parameters R ~ 4.0, (k±) 
=0.14 GeV 2 , (£ 2 ) =0.2 

Fig. 2 The same Pion Form Factor for (£ 2 ) = 0.3 
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